Initial-conditions problem for a Chiral Gross-Neveu 

system 

P.L. Natti * 

Institute de Ffsica Teorica, Universidade Estadual Paulista, 
Rua Pamplona,145 - 01405-900 Sao Paulo, S.P., Brasil 

A.F.R. de Toledo Piza 
Instituto de Ffsica, Universidade de Sao Paulo, 
CP. 66318, 05389-970 Sao Paulo, S.P., Brasil 

February 1, 2008 



Abstract 

A time-dependent projection technique is used to treat the initial- value problem for 
self-interacting fermionic fields. On the basis of the general dynamics of the fields, we 
derive formal equations of kinetic type for the set of one-body dynamical variables. A 
nonperturbative mean-field expansion can be written for these equations. We treat this 
expansion in lowest order, which corresponds to the Gaussian mean-field approximation, 
for a uniform system described by the Chiral Gross-Neveu Hamiltonian. Standard 
stationary features of the model, such as dynamical mass generation due to chiral 
symmetry breaking and a phenomenon analogous to dimensional transmutation, are 
reobtained in this context. The mean- field time evolution of non-equilibrium initial 
states is discussed. 



* Present address : Departamento de Ffsica, Universidade Estadual de Londrina, Parana, Brasil 



1 



1 Introduction 



Over the last two decades, interest in the initial-conditions problem for field theoretical 
models has been kindled and sustained by the needs of such apparently diverse 
cosmology [J] and several branches of many-body physics, notably in connection with the 
analysis of transient phenomena in the collision of complex nuclear systems 0. In these 
contexts one typically tries to obtain and solve equations describing the kinetic behavior of 
a particular, "relevant" subsystem or of a restricted set of "relevant" observables of a more 
comprehensive autonomous system. Such is the case e.g. of the scalar driving field in the 
inflationary scenario of the early universe and of one-body densities and certain correlation 
functions in heavy-ion collisions. In general, the relevant properties can be retrieved from 
appropriately constructed reduced density operators in the Schrodinger picture, which will 
evolve nonunitarily on the account of correlation effects involving different subsystems ||^ 
The nonunitary effects will manifest themselves through the dynamical evolution of the 
eigenvalues of the reduced densities, so that particular subsystems will in general evolve in a 
nonisoentropic manner [^. 

The overwhelming complexity of such a picture is considerably reduced whenever one 
is able to find physical grounds to motivate a mean-field-like approximation which consists 
in assuming isoentropic evolution of a relevant subsystem under effective, time-dependent 
Hamiltonian operators for each subsystem In this case the dynamics of the subsystem 
density matrix can be formulated in terms of a Liouville-von Neuman equation governed 
by an effective Hamiltonian and studied e.g. from the point of view of the functional field- 
theoretical Schrodinger picture, as proposed by Jackiw Unfortunately the resulting prob- 
lem still involves in general nonlinear Hamiltonians, and cannot be solved without further 
approximation. In the field-theoretical context, this has been implemented through the use of 
a Gaussian ansatz for the subsystem density functional in the framework of a time- dependent 
variational principle supplying the appropriate dynamical information, notably for bosonic 
fields §. 

It is not difficult to see that this last approximation amounts to a second mean-field 
approximation, now at the microscopic level of the single-field, nonlinear, isoentropic effective 
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dynamics. Actually, the Gaussian ansatz, having the form of a exponential of a quadratic 
form in the field operators, implies that many-point correlation functions can be factored 
in terms of two-point functions. This is well known in the context of the derivation of the 
Hartree-Fock approximation to the nonrelativistic many-body problem This factorization 
has been used by Chang |^ to implement the Gaussian approximation for the theory. The 
dynamics of the reduced two-point density becomes then itself isoentropic, since irreducible 
higher-order correlation effects are neglected. 

The focus of this work is a reevaluation of this second mean-field approximation, for 
fermionic fields, in terms of a time-dependent projection approach developed earlier for the 
nonrelativistic nuclear many-body dynamics by Nemes and de Toledo Piza p[ . This approach 
allows for the formulation of a mean-field expansion for the dynamics of the two-point correla- 
tion function from which one recovers the results of the Gaussian mean-field approximations 
in lowest order, i.e., this approach permit to include and to evaluate higher dynamical correc- 
tions effects to the simplest Gaussian mean-field approximation. Moreover, the expansion is 
energy-conserving (for closed system) to all orders The resulting dynamical equations ac- 
quire the structure of kinetic equations which eliminate the isoentropic mean-field constraint 
describing the effective dynamics of a selected set of observables P]. This approach was 
recently applied for the solution of the self-interacting A^"^ theory in (1+1) dimensions |T0[ . 



Lin and de Toledo Piza find that the Gaussian mean-field approximation fails qualitatively 
and quantitatively in the description of certain field variables. These failures are partially 
corrected by the collisional terms. Motivated by success obtained in description of time 
evolution of a off-equilibrium uniform boson (scalar field) system beyond Gaussian mean- 
field approximation in quantum-field theoretical context, it becomes interesting to study the 
fermion case in this approach. As a first step towards this end we consider in this paper 
the implementation of the Gaussian approximation to a self-interacting system of fermions. 
This is done in a framework suitable for the subsequent inclusion of collisional (correlation) 
corrections along the lines developed in ref. for the self-interacting bosonic field. 

We consider, for simplicity and definiteness, the case of an off-equilibrium, spatially uni- 
form many- fermion system described by Chiral Gross- Neveu model |jTlj|. This is an interesting 
non-trivial, renormalizable model for which many results are available in the literature so that 
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it offers suitable testing ground for the proposed approach. On the basis of the general dy- 
namics of the fields, we derive equations of kinetic type for the set of one-body variables 
in lowest order, which correspond to the Gaussian mean-field approximation. The detailed 
consideration of correlation corrections is deferred to future work. 

An outline of the paper is as follows. In Sec. II we obtain the dynamical equations which 
describe the time evolution of a general uniform fermion system. These equations are the 
groundwork for the implementation of the time- dependent projection technique. This tech- 
nique and the approximation scheme are described in Sec. III. In Sec. IV we implement in 
the quantum-field theoretical context the projection technique and obtain in Gaussian mean- 
field (isoentropic) approximation the equations which describe the effective dynamics of a 
off-equilibrium spatially uniform (1+1) dimensional self-interacting fermion system described 
by Chiral Gross-Neveu model |0. In Sec.V we use the static solution of these equations in 



order to renormalize the theory, leading to the well-known effective potential obtained by 
Gross and Neveu using the 1/A^ expansion. In this same section, we show also that other 
static results which have been discussed in the literature such as dynamical mass generation 
due to chiral symmetry breaking and a phenomenon analogous to dimensional transmutation 
can be retrieved from this formulation in a mean-field approximation. Finally we obtain 
and discuss numerical solutions for non-equilibrium initial conditions. Sec. VI is devoted to a 
final discussion and conclusions. Some points of a more technical nature are discussed in the 
Appendices. 

2 Kinetics of a self-interacting fermionic field 

In this section, we shall describe a formal treatment of the kinetics of a self-interacting 
quantum field. Although the procedure is quite general, we will adopt the specific context of 
a spatially uniform fermion system. We will illustrate all the relevant points of the approach 
and cut down inessential technical complications. Features of more general contexts are 
discussed in Ref. [Q. 

The idea of our approach is to focus on the time evolution of a set of simple observ- 
ables. We argue that a large number of relevant physical observables are one-body operators. 
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Consequently, the time evolution of observables which involve field bilinear forms such as 
iIj{x)iIj{x) , iIj{x)iIj{x) , .... is desirable. These are the observables which are kept under direct 
control when one works variationally using a Gaussian functional ansatz, and will therefore 
be refered to as Gaussian observables. In order to keep as close as possible to the formulation 
appropriate for the many-body problem, we work in fact with expressions which are bilinear 
in the creation and annihilation parts of the fields in momentum space with periodic bound- 
ary conditions in a spatial box of lenght L, defined in terms of an expansion mass parameter 
m. We begin by expanding the Dirac field operators ip{x) and in Heisenberg picture as 



where i ^k,i [&k 2 ^k,2] are fermion creation and annihilation operators associ- 
ated with positive[negative]-energy solutions 'Ui(k) [M2(k)] of Dirac's equation. Canonical 
quantization demands that the creation and annihilation operators satisfy the standard an- 
ticommutation relations at equal times 

{blx{t),bi,,,x'{t')}t=t' = (^k,k'V' for A, A' = 1,2 

(2) 

{bU{t),bl,At')}t=t' = {&k,A(i),fck',v(OW = . 



In Eq.(l) X denotes the spatial coordinate only and we use the notation 

(ko)^ — (k)^ -I- and kx — kot — kx . 

In general, the state of the system is given in terms of a many-body density operator JF 
in the Heisenberg picture, a time independent, non-negative, Hermitian operator with unit 
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trace. At this point, however, we speciahze the analysis to the case of spatially uniform sys- 
tems. These systems exhibit translational invariance (homogeneity) and rotational invariance 
(isotropy). In the case of (1+1) dimensions this reduces to invariance under translations and 
under reflection. The possible nonvanishing mean values of bilinear forms of field operators 
are 

i?k,A';k,A(t) = Tr[(6^,(t)6k,A(t))^] for A, A' =1,2 

(3) 

nk,A';k,A(t) = Tr[(6_k,A'(t)VAW)-^] for A, A' = 1,2. 



The Hermitian matrix R and the antisymmetric matrix 11 are the one-fermion density and 
pairing density respectively. Using these objects we can construct the extended one-body 
density [|T^ 



i?k(t) nk(t) 
-nut) h-Riit) 



L {bUyit)bUt)) (&k,A'(t)6LW) 



(4) 



This object summarizes all information on the Gaussian observables and provides an adequate 
starting point for our kinetic treatment. 

The first step is standard and consists in reducing the extended one-body density to 
diagonal form. This can be achieved by subjecting the creation and annihilation operators to 
a canonical transformation of the Bogolyubov type with coefficients defined by the eigenvalue 
problem 



Af^(^)7^k(t)A'k(t) = Qk{t) 



(5) 



where the unitary matrix X^t) which diagonalizes 7^k(^) has the structure 



Y\^ Xu 



^k 



^k" -^k 



(6) 
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and Q]i{t) is a diagonal matrix which can be written as 



h- Mt) 



(7) 



The unitary conditions for A'k(t) can be interpreted as orthogonahty and completeness re- 
lations for a set of natural orbitals which diagonalize the extended one-body density. They 
read 



A'j^A'k — /4 and AkA"]^ — I4 , 
or, more explicitly, in terms of the submatrices and Yk [see Eq.(6)] 



(8) 



(9) 



Y^Y]^ + X'^X^ — I2 



Y^X* + XlY^ = O2 



The elements of the diagonal submatrix ^'k(i) can, on the other hand, be interpreted as 
quasi-fermion occupation numbers for the paired natural orbitals. Because of the assumed 
symmetry one must have 



(10) 



Finally, using Eq.(6), we can relate the fermion operators blA't) ^k,A(i) to the new 
quasi-fermion operators Ptxi^) Py^A^) for A = 1,2 as 







' Xu 


X21 




^21 






bk,2 




X12 


X22 


•y* 


-'22 




/9k,2 






Yn 


Y21 


X*n 


^21 






b\2 . 




. ^12 


Y22 


^12 


-^22 . 




. /3lk,2 . 



fill 



With the help of Eq.(ll) it is then an easy task to express and i^ix), Eq.(l), in term 
of Pl^xi^) ^^'^ Pk,x{t) foi' ^ = IjS. In doing so, one finds that the plane waves of 
and iIj{x) are modified by a complex, momentum-dependent redefinition of m involving the 
Bogolyubov parameters. The complex character of these parameters is actually crucial in 
dynamical situations, where the imaginary parts will allow for the description of time-odd 
(i.e., velocity-like) properties. 

What we have achieved so far amounts to obtaining an expansion of the fields and 
ip{x) such that the mean values in of Gaussian observables are parametrized in terms of the 
natural orbitals Xx',x{yi} and Yx'^x{\s.) and of the occupation numbers ^'k,A(^) — Tr (/3k,A/^k,A-^) 
for A = 1, 2. In general, all these quantities are time dependent under the Heisenberg 
dynamics of the field operators, and we now proceed to obtain the corresponding equations 
of motion. Taking the time derivative of Eq.(5) and using the unitarity condition (8) we get 



(12) 



We next evaluate the left-hand side of this equation using the Heisenberg equation of motion 
to obtain 



Tr ( [(3ly(3^.,x, h] t) Tr ( [/?_k,A'/?k,A, H] T) 
Tr ( [/3lk,A'4,A, h\ t) Tr ( [/3k,A'4,„ h\ T) 



(13) 



The right-hand side of Eq.(12) can also be evaluated explicitly using Eqs.(6) and (7): 



i (Qk — '^k'^kQk — Qk'^k'^k) — 



iz>k + K, /ik]- -9i + {'^k, ^k}+ 

-t/k + {i^k, 5'k}+ -i2>k + [i^kj/ik]- 



(14) 



where the matrices /ik and g^k are given in terms of Xk and Ik as 



/ik = -i{Y^YC + XlX^) 



(15) 
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From Eqs.(13) and (14) we obtain dynamical equations which describe the time evolution of 
our uniform fermion system. They read 

(16) 

Eqs.(16), together with the unitarity conditions (9), determine the time rate of change of 
the Gaussian observables in terms of expectation values of appropriate commutators. They 
are, however, clearly not closed equations when the Hamiltonian H involves self-interacting 
fields. In this case, in fact, the time derivatives of the Gaussian observables are given in 
terms of traces which are not expressible in terms of the Gaussian observables themselves, 
since they will involve also many-fermion densities. This situation can be dealt with in terms 
of the projection technique reviewed in the next section. 

3 Projection technique and approximation scheme 

In this section we introduce the time- dependent projection technique |^ which permits to 
obtain closed approximations to the equations of motion (16). It has been developed earlier 
in the context of nonrelativistic nuclear many-body dynamics and was recently applied in 
the quantum-field theoretical context to the self-interacting theory in (1+1) dimensions 
p!0| . It allows for the formulation of a mean-field expansion for the dynamics of the two-point 
correlation function from which one recovers the results of the Gaussian mean-field approxi- 
mations in lowest order. If carried to higher orders it allows for the inclusion and evaluation 
of higher dynamical correlation corrections to the simplest mean-field approximation. 

In the specific context of the equations of motion (16) we begin by decomposing the full 
density JF as 
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r = :F,{t) + r{t) , 



(17) 



where J^o{t) is a Gaussian ansatz which achieves a Hartree-Fock factorization of traces in- 
volving more than two field operators. The Gaussian density JFo(t) is chosen as having the 
form of a exponential of a bilinear, Hermitian expression in the fields normalized to unit trace 
0. In the momentum basis, it reads 



^0 



exp 


E(ki,k2) ^ki,k2&ki^k2 + ^ki.ka^ki^L + C'ki,k2&ki&k2 




Tr |exp 


I](ki,k2) ^ki,k2^ki^k2 + -Ski,k2^ki^k2 + C'ki,k2^ki&k2 


} 



The parameters in Eq.(18) are fixed by requiring that mean values in JFq of expressions that 
are bilinear in the fields reproduce the corresponding JF averages [see Eqs.(20) below]. JFq is 
a time-dependent object, which acquires a particularly simple form when expressed in terms 
of the Bogolyubov quasi-fermion operators 



^o(t) = n [^k,A4,A(t)/5k,A(t) + (1 - ^k,A)/9k,A(t)4,A(^) 



(19) 



k,A 



This is clearly a unit trace object which in addition satisfies 



Tr (/3,^o) = Tr (/5„^) = Tr (/^t^o) = Tr (/Jt^) = ; 
Tr(/3,/5fe^o) =Tr(/5,/3fe^) =0 ; 

Tr {M:F,) = Tr (/Jt/^t^) = ; (20) 
Tr {(3i(3tJ'o) = Tr = iyJa,t and 

Tr {mI^o) = Tr (mI^) = (1 - '^a)Sa,b • 



Correspondingly, the "remainder" density J^'{t), defined by Eq.(17), is a traceless, pure 
correlation density. As already remarked, a crucial point to observe is that J^oit) can be 
written ClS db time- dependent projection of J-", i.e.. 
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jCo(t) = ■p{t)r with v{t)v{t) = v{t) 



(21) 



In order to completely define this projector we require further that it satisfies 

z^o(t) = [P{t), li\T= [^o(t), ii\ + ViS) [H, ^] , (22) 
where C is the Liouvillian defined as 

C- = [H,-], (23) 

H being the Hamiltonian of the field. Eq.(22) is just the Heisenberg picture counterpart 
of the condition V{t)J-' = which has been used to define V{t) in the Schrodinger picture 
0. It is possible to show that conditions (21) and (22) make V{t) unique and to obtain an 
explicit form for this object in terms of the quasi-fermion operators and of the natural orbital 
occupations 0, |T^. 

The existence of the projector V{t) allows one to obtain an equation relating the corre- 
lation part J^'(t) to the Gaussian part J^o(t) of the full density. This can be immediately 
obtained from Eqs.(17), (21) and (22) and reads 

{zdt - V{t)C) J'' it) = (X - Vit)) CJ^oit) . (24) 
This equation has the formal solution 

J^'it) = Git, O)J^'(O) - I t dt' g{t, t') (J - V{t')) CToit') , (25) 

Jo 

where the first term accounts for initial correlations possibly contained in JF. The object 
g{t,t') is the time-ordered Green's function 



git,t')=T(exp 



I C dTV[T)C 



(26) 
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We see thus that J-''{t), and therefore also JF [see Eq.(17)], can be formally expressed in 
terms of J-'o{t') (for t' < t) and of initial correlations JF'(O). This allows us to express also the 
dynamical equations (16) as functionals of J-'o{t') and of the initial correlations. Since, on 
the other hand, the reduced density J^o{t') is expressed in terms of the one-fermion densities 
alone, we see that the resulting equations are now essentially closed equations. Note, however, 
that the complicated time dependence of the field operators is explicitly probed through the 
memory effects present in the expression (25) for J^'{t). Approximations are therefore needed 
for the actual evaluation of this object. A systematic expansion scheme for the memory 
effects has been discussed in Refs. |, |l^. The lowest order correlation corrections to the 
pure mean field approximation, in which JF' is simply ignored, correspond to replacing the 
full Heisenberg time-evolution of operators occurring in the collision integrals by a mean-field 
evolution governed by 

Ho = V\t)H . 



Consistently with this approximation, £ is replaced in (25) and (26) hj Cq • = [Hq, ■ ]. In 
this way correlation effects are treated to second order in H in the resulting collision integrals. 

An important feature of this scheme (which holds also for higher orders of the expansion 
0) is that the mean energy is conserved, namely 

where 

{H) = Tr HJ^oit) + Tr HJ^'{t) . 

In the following sections we apply the general expressions obtained in above to treat a uni- 
form fermion system described by Chiral Gross- Neveu model (CGNM). We will consider only 
the lowest (mean-field) approximation, corresponding to T'{t) = 0. Collisional correlations 
will be treated elsewhere. 
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4 The Chiral Gross-Neveu model (CGNM) 



The Hamiltonian density for the CGNM is given by 



A'' 2 



i=l 



N 



.i=l 



TV 



E^^75^^ 



i=l 



(27) 



where ^ is a constant which indicates whether the model is invariant under the discrete 75 
transformation = 0) or under the Abehan chiral U{1) group = 1). 

This is a massless fermion theory in (1+1) dimensions with quartic interaction. The 
model contains N species of fermions coupled symmetrically, where ip^ is a complex Dirac 
spinor transforming as the fundamental representation of SU (N) group. It is known that the 
actual symmetry of the theory is not SU{N) but rather 0{2N) [|l^. The transformations 
forming this group mix not only particles but also particles with antiparticles. This model 
is essentially equivalent to the Nambu-Jona-Lasinio model [|15|, except for the fact that in 
(1+1) dimensions it is renormalizable. Moreover, it is one of the very few known field theories 
which are assimptotically free. To leading order in a 1/A^ expansion |]TT|, the CGNM exhibits 
a number of interesting phenomena, like spontaneous symmetry breaking [l^, dynamical 
fermion mass generation and dimensional transmutation. The model possess an infinite 
number of conservation laws, and as a consequence, the S'-matrix may be computed exactly 



To obtain the time evolution of Bogolyubov parameters we have to obtain the CGNM 
Hamiltonian [see Eq.(16)] by integrating over one-dimensional space. This involves, in par- 
ticular, choosing a representation for the 7-matrices. Here we have to be careful, since a bad 
choice of representation can spoil manifest reflection invariance. In Appendix A (see also 



Ref. ||15||) we discuss suitable representations for the 7-matrices. We choose the Pauli-Dirac 
representation, namely 



7o = 0-3 



7i = ia2 and 75 = 7071 = ai 



(2J 



In this representation the spinors ui(k) and U2(k) are given by 
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(fcp + m) 
2m 



1/2 



1 

k 



. (ko + m) _ 



, 1*2 (k) 



{kg - m) 
2m 



1/2 



1 

k 



(/co - rn) 



(29) 



The resulting form for the CGNM Hamiltonian obtained using this representation is given in 
full in Appendix B. 

We next consider the initial conditions problem for this system. The general Bogolyubov 

transformation defined in (11) breaks both chiral and charge symmetries, but we restrict the 
following development to a special Bogolyubov transformation (to be called Nambu transfor- 
mation) which breaks the chiral symmetry of our system only. The elements of this Nambu 
transformation, parametrized consistently with the unitarity conditions (9), are given by 

Xii = X22 = cos <^k and X12 = ^21 = 

(30) 

Yi2^-Y2i^sm(fi^e'''^ and = I22 = . 



In the special case of a Nambu transformation, the elements of the matrices and [see 
Eq.(15)] are given by 



^11 = ^22 = 7k sin^ 

^12 = /i21 = S'll = 5'22 = (31) 

5'i2 = -5'2i = [i(Pk - 7k sin (/?k cos (/7k] e*'^'' . 



On the other hand, in the mean- field approximation one finds that ^o(^) commutes with 
number operators 

Tr {i/cGNM[J^o, PliM} = Tr {Hca^^J^o, /3I2/3K2]} = , (32) 
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while, due to the charge conservation in the system, we have also 



(33) 

Tr {[/?-k,i/3k,i, -f^cGNMl^o} = Tr {[/?-k,2/3k,2, -f^coNMl^o} = . 



Substituting Eqs. (31), (32) and (33) in Eq. (16), we obtain the equations which describe 
the time evolution of our system as 

z>k,i = and z>k,2 = (34) 
bV'k + 7k sm (^k cos (pk\ e = — ^ '■ r . (35) 

U — TAi,l — i^k,2j 



Equation (34) shows that the occupation numbers of the paired natural orbitals are con- 
stant, i.e., we recover the general isoentropic character of the mean-field approximation. The 
complex equation of motion (35) describes the time evolution of the Nambu parameters. 
Writing the CGNM Hamiltonian, given in Appendix B, in the Nambu basis using Eqs. (11) 
and (30), and substituting this Hamiltonian in Eq.(35), we obtain the explicit form of these 
equation. The calculation of traces is lengthy but straightforward. As a result one obtains 



. sin 2</?k (k 



+ 



2 

(e + 1 
47r ) ko 



sin 2(^k — m— 
ko ko 



sin^ (^kC ''^^ - cos^ ^Pkc''^^ 



sm . 



'^'Pk + — f sin^ (p^e '^"^ - cos^ (Pkc'^'' 
m ^ 



+ 



(/i + h 



(36) 



where Ii and I2 are the divergent integrals 



15 



r dk' 

= / — cos 2</7k'(l - - ^'k',2) 

J Kq 



f dk! Ik' I 

h = -r, sin 2(^k' cos 7k' (1 - f^k',i - • 



(37) 



We take = 1 for simplicity. Splitting the complex equation (36) into real and imaginary 
parts we have 



(/7k = sm7k 



7k sin 2(/3k 



2 sin 2(/3k 
ko 



47r 



+ 



+ 2 cos 2(/7k cos 7k 



ko 



(38) 



Finally, the mean-field energy is evaluated as 



/ ffM.F. \ 
V-'-'CGNM/ 



+ 



(39) 



87r2 ; 2 



(ik'(i + i/k',1 - i^k',2) <ik"(i - i/k",i + i^k",2) , 



where Ii and /2 are given in (37) and is given by 

f dk' /k'\^ 
^3 = / TT ~ cos2(/3k(l - - I^k',2) • 



(40) 



Since all the results above contain divergent integrals, a renormalization procedure is 
required. This is discussed in the next section. 
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5 Renormalization 



In order to handle the infinities which occur in the preceding equations it is necessary to 
introduce a renormahzation prescription that will render physical quantities finite. In gen- 
eral, renormalization procedures consist in combining divergent terms with the bare mass 
and coupling constants of the theory to define new, finite (or renormalized) values of these 
quantities. Thus, when one adopts a momentum cut-off to regularize divergent integrals, the 
bare mass and coupling constants are chosen to be cut-off dependent in a way that will cancel 
the divergences. In the present case, however the divergent integrals (37) and (40) involve 
the dynamical variables themselves in the integrand, so that even their degree of divergence 
is not directly computable. In order to handle this situation we will use a self-consistent 



renormalization procedure inspired in Ref.[18 



The renormalization prescription we use can be based on the consideration of the static 
solutions of the dynamical equations (38). They are determined by the equations 



sm7k|cq 



tan 2y9k 



eq 



-|k|m[l-(^V47r) (^ + !)(/, + /,)] 
[(k)2 -f (^W/47r) (e + l)(/i + /2)] 



cos 7k I eq • 



(41) 
(42) 



We will then show explicitly that it also controls the divergences which appear in the kinetic 
regime of the mean-field approximation. 

In order to obtain the renormalization prescription we introduce a momentum cut-off A 
and begin by assuming that, in order to render the theory finite, the bare coupling constant 
must approach zero for large values of A as (see e.g. Refs. []rT|, |T^) 



Ait 



K + i) 



n -1 



(43) 



where the form of the first factor is dictated by later convenience. We next assume that the 
integrals Ii and I2 have logarithmic divergences 
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h = a + Mnf^) (44) 



h = c + dlnf^l , (45) 



where a, 6, c and d are finite constants. Substituting (43) and the ansatze (44) and (45) in 
the static equation (42) we obtain 

tan 2^ I _ -(-irHk|[l-(& + rf)] 

- [(k)2 + m^(6 + d)] ■ 

Finally, evaluating the integrals /i and /2, Eqs.(37), using Eq.(46), we are able to determine 
the constants a, b, c and d self-consistently. This calculation is given in Appendix C. We find 
there b — 1 while d remains arbitrary. The renormalized static equations are then obtained 
by simply substituting these values in (46) . We obtain 

„ , (-l)"m|k|d 

We observe that the theory involves just one free parameter, d. This is altogether reasonable 
since our starting point was a massless fermion theory which was determined by one dimen- 
sionless coupling constant g"^, and we end up with a theory determined by one free parameter 
d after the self-consistent renormalization procedure. What is needed is an interpretation of 
the parameter d. 

To obtain this, we begin by writing the fields ■0(x) and ■0(x) given in (1) in the Nambu 
quasi-particle basis using (11), and (30). We find that the new Dirac spinors in this basis are 
given by 



u[(k.) — cos</7k Ui(k) + sin (/7k e^'^ (—k) 

(48) 

M2(k) = cos(pk 'W2(k) — sin e~'^'^^ mi(— k) . 



18 



Prom the renormalized static solution (47) we have that 



COS(/?k|eq = 



Sm(/?k|eq = 



1 



{kox + m|k|d)^/2 + {kox - m\k\dy/'^' 



(49) 



{koX + m|k|c?)^/^ - {koX - m\k\dy 



/2' 



where x — [k'^ + {1 + d)^m^]^/^. Substituting the solutions (49) and the particle spinors (29) 
in (48), we obtain the new spinors in quasi-particle basis. They are given as 



^i(k) = 



«2(k) 



.1/2 



2m. 



eff 



2m, 



eff 



1/2 




(50) 



where 



{ko^-'^y = {kf + {meffY and me// = (1 + d)m 



(51) 



Comparing the spinors (50) in quasi-particle basis with the spinors (29) in particle basis, we 
see in Eq.(51) what amounts to a redefinition of the mass scale. Therefore, we note that, 
unlike the situation found in connection with the 1/N expansion, the use of the Gaussian 
ansatz, Eq.(19), parametrized by the canonical transformation leading to the quasi- fermion 

basis, allows for the direct dynamical determination of the stable equilibrium situation of the 
system [see Eqs.(41) and (42)], including symmetry breaking and mass generation. Moreover, 
the renormalization procedure effectively replaces the dimensionless coupling constant 
by the free parameter d associated to the mass scale [see Eq.(51)]. This is analogous to 
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the phenomenon of dimensional transmutation found by Gross and Neveu [|TT| in the 
expansion. Finally, aside from the over-all mass scale (characterized by d) there are no free 
adjustable parameters. 

Using Eqs.(38) and (43)-(45), we finally write the renormalized form of the dynamical 
equations that describe the mean-field time evolution of this system. As mentioned before 
they are now also finite and read 



7k sin 2(pk 



-md— sin 7k 



(52) 



2 sin 2(y9k 
kn 



{\if + m\l + d) 



+ 



— 2md— cos 2y3k cos 7k . 

The first two equations express just the isoentropic character of the mean- field kinetics. The 
remaining coupled equations describe the kinetic behavior of the Nambu parameters and are 
discussed in the following section. 

We conclude this section calculating the ground-state (vacuum) energy of our system in 
the mean-field approximation. This approximation to the ground state is obtained as a static 
solution of Eqs.(52). The associated mean energy can be obtained by taking z/k,i = z/k,2 = 
in (39) and evaluating the divergent integrals (37) and (40) using the renormalization ansatze 
(43) and (44)-(45). We obtain 



/ ZTM.F. \ vacuum 
\ CGNM/ 



1 

2^ 



1 

2^ 



:i + d)v 



In 



[l + df 



(53) 
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where we used for the finite constants a, b and c the expresions given in Appendix C. The 
first term, divergent as A ^ oo, represents a vacuum background energy. The presence of 
this term has no physical consequences, and we follow the usual practice of redefining the 
zero of the energy scale by simply subtracting it out. Therefore, in terms of the effective mass 
given by (51), the renormalized energy density of the ground-state of our system is given by 



/ ZJM.F. \ vacuum 
n-'^cgnm/ 



m 
An 



m. 



m 



l + 21n2-ln 



eff 



m 



(54) 



Fig.l shows the renormalized ground state energy density as a function of effective mass 
rrieff. This figure reproduces the well-know effective potential obtained in the case of 1/A^ 
expansion [|ni|. 



6 Kinetics of One-Fermion Densities 

In this section, we discuss the solutions of the nonlinear equations of motion of the Nambu 
parameters and 7k in terms of the values taken by the parameter d. In order to do this, 
it is useful to note that these equations can be obtained as the canonical equations of the 
effective c-number Hamiltonian 



hcs 



2m|k| 
kn 



(icOS7k(l — CTy^ 



2^1/2 



(55) 



where we defined = — cos2(y9k as a new momentum-like variable canonically conjugate 
to 7k. For a general value of d this dynamical system has equilibrium solutions given by 
nn with v^kjeq given by Eq. (47). A special situation occurs, however, when d = 0. 



7k 



eq 



In this case, in fact, y^k becomes a constant of motion (7k becomes cyclical in h^g) and 7k 
becomes independent of ip^. 

Considering next the time-dependent solutions in the general case {d ^ 0), we see that 
the equilibrium solutions are stable in the sense that small initial displacements from the 
equilibrium states lead to librational behavior both in y^k and in 7k. For sufficiently large 
displacements, however, the latter variable aquires a rotational (though nonuniform) behavior 
while v^k (or c^k) still oscillates about an equilibrium value. A set of phase-space trajectories 
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illustrating this behavior is shown in Figs. 2 and 3. As one sees in these figures, when d 
is closer to the special value zero, the stability domains shrink in the variable as the 
equilibrium points approach the lines (/^^ = m7i/2. If, on the other hand, one considers the 
kinetic equations for d = 0, one sees that the solutions correspond to uniform 7k rotations 
with fixed <^k- 

In order to interpret this behavior we recall that d, together with the expansion mass m, 
and the equilibrium values of (pk and 7k, define the equilibrium effective mass me// [see Eq. 
(51)]. In particular, when d — —1, one has rrieff — and the equations of motion reduce to 



Ikl 

0k = m— sm7k 



7ksin29?k = 2— — sin2<^k + 2m— cos2v9kCOS7k 
ko ko 



which reproduce the equations one obtains for a free (massless) field e.g. by making g — 
in the Hamiltonian density (27). Other values of d will, on the other hand, correspond to 
nonvanishing equilibrium values of the effective mass. Displacements of (/?k and 7k from their 
respective equilibrium values will then correspond to preparing the system in a state having 
a "wrong" (k-depcndent) effective mass (in particular, non-vanishing in the "symmetric" 
phase, d = —1). The ensuing motion represents therefore the dynamical reaction of the 
system to this situation. Several phase space trajectories for off-equilibrium solutions are 
shown in Figs. 2 and 3. 

The special case d — coresponds to choosing the equilibrium effective mass mg// itself 
as the expansion mass m. In this case the stability islands disappear, as the phase space tra- 
jectories degenerate to straight lines. The off-equilibrium dynamics can then be represented 
by the time dependence of the single parameter 7k. 
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7 Discussion and conclusions 



We have described a treatment of the initial- value problem in a quantum field theory of self- 
interacting fermions in the Gaussian approximation . Although the procedure is quite general, 
we have implemented it for the vacuum of a relativistic many-fermion system described by 
Chiral Gross-Neveu model (CGNM). 

We have obtained the renormalized kinetic equations which describe the effective dynam- 
ics of Gaussian observables in the mean-field approximation for a uniform (1+1) dimensional 
system. We used the static solution of these equations in order to renormalize the theory, 
leading to an effective potential similar to that obtained by Gross and Neveu using the 1/A^ 
expansion. We show also that other static results discussed in the literature such as dy- 
namical mass generation due to chiral symmetry breaking and a phenomenon analogous to 
dimensional transmutation can be retrieved from this formulation in the mean-field approxi- 
mation. Finally, we obtained and discussed numerical solutions for the time evolution of the 
Nambu parameters (Gaussian variables) for non-equilibrium initial conditions. 

As a final comment we note that, unlike the situation found in connection with the 1/A^ 
expansion, the use of the Gaussian ansatz, Eq.(19), parametrized by the canonical trans- 
formation leading to the quasi-fermion basis, allows for the direct dynamical determination 
of the stable equilibrium situations of the system [see Eqs.(41) and (42)], with or without 
symmetry breaking and mass generation. In the latter case we obtain just a free, massless 
theory. 

The present projection technique can be extended to include and to evaluate higher 
dynamical correlation corrections to the mean-field approximation. In this case the occu- 
pation numbers are no longer constant, z>k,A 7^ 0, and their time dependence affects the 
effective dynamics of the Gaussian variables (see reference [0). A finite matter density 
calculation beyond the mean-field approximation allows one to study collisional observables 
such as transport coefficients [^. Finally, we comment on the extension to non-uniform 
systems. In this case the spatial dependence of the fields ipi^x) and ipi^) are expanded in 
the natural orbitals through the use of a non-homogeneous Bogolyubov transformation (see 
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reference [9]). 



Acknowledgments 

One of the authors (P.L.N.) was supported by Conselho Nacional de Desenvolvimento 
Cientifico e Tecnologico (CNPq), Brazil; and by Fundagao de Amparo a Pesquisa do Estado 
de Sao Paulo (FAPESP), Brazil. 

Appendix A : Uniform system and the representations 

for the 7^- mat rices 

A (1 + 1) dimensional uniform system has to be invariant under translations and under re- 
flection. Reflection invariance implies that the equations of motion are direction independent, 
i.e., dynamical quantities should involve |k| only. 

Let us consider the solutions of the Dirac equations 



(i^-m,-)V'^-^H^) = 0, J = 1,2 , 



(56) 



for two different values mi and m2 of the mass: 



mi 



1/2 



^('"i)(x) = 



nil 



1/2 , 



(57) 

(58) 



and 



m2 



1/2 



(m2) 



(59) 
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7712 



1/2 



(m2) 



(60) 



By imposing (say at t — 0) the condition 



(61) 



one can show that the operator sets (6k|i^\ b^2^) and {b^i\ b^2^) are related by the canonical 



transformation 



(™2) 



1712/ 



("12 ) 





1/2 



«i'"^nk)wj™^^(k) 



+ < 



1712/ 



(m2) \ 




/in 



xzi"^^^(k)iif^\-k) 



^6 



(mi)t 
-k,2 



(62) 



^k,2 



— < 



1/2 



(mi) 



m2/ 



1/2 / t 



("12) " 



1/2 



tzr^)(k)x.^)(k) 



(-k)«r^(k) 



^-k,l 



(63) 



Reflection invariance requires that the coefficients of this canonical transformation be 
invariant under k — > — k, i.e. 



u^r\^)ut^'\^) = «r'(-k)wr^(-k) 



71 (■'■"2)/ 
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(64) 



There are just two representations of the 7'^-matrices satisfying these conditions. These 
are the Pauh-Dirac representation 



To = 0-3 ; 7i = and 75 = 7071 = (7i 



(65) 



and 



lo = CTs ; 7i = *cri and 75 = 7071 = -a2 . 



(66) 



We choose the Pauh-Dirac representation (65) which admits as solutions the spinors given in 
(29). It should be noted, however, that the spinors 



2m 



1/2 



1 
k 



(ko + m) 



U2ik) 



' ko + 
2m 



1/2 



{ko + m) 
1 



(67) 



which are also solutions of Dirac's equation in the Pauli-Dirac representation, do not satisfy 
the uniformity conditions (64). 



Appendix B : The uniform CGNM Hamiltonian 

In this appendix we give explicitly the Hamiltonian of the Chiral Gross-Neveu model in a 
form which preserves the reflection invariance of the system. In the Pauli-Dirac representation 
(see Appendix A) the CGNM Hamiltonian is given in full by 
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N ^ 



(9_^)y y I 





^k-M^k'" 




1 [esi(k075"i(kn^ii(k;")75"i(kf ) 




ai(k;")7ii(kf )] 


'^k^+ky.k^'+ky" 








[eui(k;)75^ii(kn^x2(k;" 


)75'U2(k"") 


- ■ui(k^?ii(k'/ 


\ - /] /// 


)U2{kj ) 


^k;+k^"',k;'+ky' 




i^k'.',i 
» "-J 


,i^k"' 


[^«i(k075«i(kn^i(k7 


)75«2(kf) 


-wi(k^)txi(k7 


)MK 


)«2(kf)] 


* J 3 ^ 


+ bl 


i^k'.'.iV." 


,2^k"' 


[^ui(k075«i(kn^i2(kf 


)75«i(kf) 


-nl(k0ul(kf 


)MK 




* * J 3 


+ 6k; 


2^k;',24^" 




[^U2(k075W2(knwi(kf 


)75^^i(kf) 


- M^M^i 


)MK 




(^k'+k'/",k'/+k'" 

t 3 ' I ' 3 


+ V. 






, 2 [^n2(k075t^2(kn«2(kf 


)75t^2(kf) 




)MK 




^k'+k'".k"^k"" 

^ J 3 


+ 


2&I",2felc"' 


A'" 

J 


, 2 [C«2(k075«2(kntii(k;" 


)75«2(kf) 


- M2(k;)u2(kf 






^k',k"+k"'+k"" 


+ fok; 


2&L',2^k;" 


.2&k"' 

J 


M [eS2(k^)75W2(knS2(k;" 


)75«l(kf) 


-U2(kDw2(k:' 


)U2{K 




'5k'.+k"'+k"",k" 




l&I",2&lc"' 






)75«l(kf) 


-^Il(k0^i2(k'/ 


)si(k;" 


)«l(kf)] 


'^k'+k'/.k'/'+k"" 


+ 6t. 




A'" 


, , [^Mk',)j,U2{K)u2iK 


)75«2(kf) 


-Ul(kO,i2(k'/ 


)^2(k;" 


)«2(kf)] 


'5k'.+k"+k'."',k"' 


+ bl 


lV.',2°k'." 


' 3 


[^«i(k;)75«2(k;')«i(k7 


)75«2(kf) 


-«i(k;)tt2(k7 


)MK 




^k^+k^'+ky+kv" 


+ 6t. 


i^k;',2^" 




[^Ui(k075W2(knW2(kf 


)75ni(kf) 


-ni(k:)u2(k^' 






<^k^+k;',ky'+ky" 


+ 6k^ 


2&k'.',l^'i,/// 


At 


[^U2(k075wi(kn«i(kf 


)75t^i(kf) 


-n2(kO«l(kf 


)MK 




^k'+k"+k'."',k"' 


+ V, 


2&k'.',lV" 


Ak'" 

3 


[^n2(k075^^l(kr)n2(k;" 


)75^^2(kf) 


-n2(k;)ui(kf 


)MK 


)MK')] 


<5k'+k"+k"',k"" 


+ &k; 


2&k",l&k'" 


.A'" 

3 


, 2 [?S2(k075«i(kn«i(kf 


)75^^2(kf') 


-«2(kOMi(kf 






^k'.+k",k"'+k"" 


+ &k^ 


2^k'/,l^k'." 


,2^k'/' 

' 3 


.,1 [ew2(k^)75«i(kn«2(k;" 


)75Wi(kf) 


-U2(kOwi(k^' 


)U2{K 




^k^+kj'+ky+k^'" 



This Hamiltonian is written in a particle basis. In the dynamical equations (35) we 
have to use the CGNM Hamiltonian in the Nambu quasi-particle basis. The transformed 
Hamiltonian is obtained in a straightforward way using (11) and (30). 



Appendix C : Self-consistent renormalization 
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Using Eq. (46) we obtain sin2</?k|cq and cos29?k|eq as 



Sin2(/7k|eq 



COs2(^k|eq = 



-(-l)"m|k|[l- (6 + ci)] 
/to[A;2 + m2(6 + d)2]V2 

[A;2 + m2(6 + d)] 
A;o[A;2 + m2(6 + d)2]V2- 



(69) 



Substituting (69) in 7i and /2 given in (37), taking i/k^i = i/k,2 = 0, and performing the 
integration, we have 



h = 



dk [k2+m2(6 + d)] 

= a + om 



(P + m^) [fc2 +TO2(6 + d)2] 



Vm2 y 



2d 



[(l + rf)2-l]V2 

when (1 + rf)^ > 1 



arctan[(l+d)2-l]i/2 + ln 



(l + rf)2 



and 6=1 



a 



1 



[l-(l+d)2]l/2 

when (1 + d)^ < 1 



1 + (l + rf)2}l/2 



l-{l-(l + (i)2}l/2 



In 



(l + d)2 



and 6=1 



(70) 



(6 + d - 1) J 



dk 



A (P + m2) [fc2 + m2(6 + d)2]V2 



c + din 



/A2X 
Vm2y 



-2d 



[(l+d)2-l]l/2 

when (1 + d)^ > 1 



arctan[(l + d)2-l]i/2 + (;in 



(l + d)2 



and d = d 



-d 



In 



[1- (l + d)2]V2 

when (1 + d)^ < 1 



l + {l-(l + d)2}l/2 



1-{1- (l + d)2}l/2 



+ din 



(l + d)2 



and d = d 



(71) 



This shows that the ansatz (44)-(45) is self-consistent. Eqs. (70) and (71) give the values 
of the constants a, h and c. We observe that d remains arbitrary in this calculation. 
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Figure Captions 



Figure 1 - Mean-field effective potential for the ground-state of an uniform fermion system 
described by the CGNM as function of effective mass rueff- 

Figure 2 - Phase-space of Nambu parameters for |k| — m — y/2 in the case of free fermion 
system [{d+1 — 0) or me// = 0]. Initial- values: 

Dotted line : ip"^ = (-0.2 + n7r/2) , -f^ = 2n7r and ip^ = (+0.2 + n7r/2) , Yk = (2^ + l)7r, 
Dashed line : ip'^ = (-0.1 + n7r/2) , Yk = 2n7r and ip'^ = (+0.1 + n7r/2) , Yk = (2ri + 
Solid line : (^jf = (-7^/4 + n7r/2) , 7if = and <^i,^ = (0 + mi/2) , 7jf = 0, 
Dot-Dashed line : ip'^ = (tt/S + nn/2) , -f'^ = 0. 

Figure 3 - Same as Fig. 2 with {d + 1) = 0.5858. Initial values: 
Dotted line : <^j,° = (-0.1 + n7i/2) , = 2mi and v^j," = (+0.1 + mi/2) , -f'^' = (2n + l)7r, 
Dashed line : (^jf = (-0.05 + n7r/2) , Y^ = 2n7r and p}'^ = (+0.05 + n7r/2) , Y^ = (2n+ l)7r, 
Solid hne : = (-tt/S + n7r/2) , tL"" = and (^jf = (0 + n7r/2) , 7]^ = 0, 
Dot-Dashed line : pi';^ = (37r/16 + n7r/2) , 7";^ = 0. 
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